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Abstract 

A new model is presented for simulating coherent synchrotron radiation (CSR) in one dimension. 
The method is based on convolving an integrated Green function (IGF) with the longitudinal 
charge density. Since it is based on an IGF, the accuracy of this approach is determined by how 
well one resolves the charge density and not by resolving the single particle wake function. Since 
short-range wakcficld effects arc included analytically, the approach can be much more efficient 
than ordinary (non-IGF) approaches in situations where the wake function and charge density 
have disparate spatial scales. Two cases are presented: one derived from the full wake including 
short-range effects, and one derived from the asymptotic wake. In the latter case the algorithm 
contains the same physics as others based on the asymptotic approximation, but requires only the 
line charge density and not its derivative. Examples are presented that illustrate the limitations of 
the asymptotic-wake approximation, and that illustrate how microbunching can increase the CSR 
fields by orders of magnitude depending on the microbunching wavelength. 



I. INTRODUCTION 

Coherent synchrotron radiation (CSR) is one of the most important and difficult to model 
phenomena affecting future light sources and other lepton accelerators involving high inten- 
sity, high energy, short bunches. It is also an important phenomena in certain astrophysical 
settings. Analyses of the radiation from a charged particle moving on a circle of radius p 
were performed by many authors including Schott [IJ and later by Goldreich and Keely [2] 
and by Murphy, Krinsky, and Gluckstern [3]. As shown in [3], at relativistic velocities, in 
steady state, in the small angle approximation, and assuming free space boundary condi- 
tions, the single particle wake function corresponding to the radiation component of the 
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azimuthal electric field is, in MKS units, w = j— where the normalized wake 



function is. 



and where 



w{p) = (1) 



In the above, Vl = p + ^J p? + 1, P = '^s, and s is the arc length between the observa- 
tion point and the source charge. Note that this is the wake for /i > 0; in the large 7 
approximation the wake is zero for yU < 0. Note that for small argument, w has the form, 

w = 1 — y/i^H , and, for large argument, w ~ — + f/U^^H • The single particle 

wake, w, is plotted as a function of s in Fig. [T] along with its asymptotic approximation. 
Given the asymptotic form of w and the fact that /i oc 7'^ it is clear from these formulas 
(which are based on the large-7 approximation) that analyses involving just the asymptotic 
wake produce CSR fields that are independent of 7. In contrast, the short-range wake does 
produce energy-dependent effects. 

The simulation of one-dimensional CSR effects usually involves computing the convolution 
of a wake function, w, with a numerically computed line charge density. A, or its derivative 
dX/ds. If the single particle wake function depends only on the separation between the source 
and observation point, then the radiation electric field can be expressed as the convolution, 

E{s) = [ dz'\{s')w{s- s'), (3) 



or 



Eis) = - I ds'^uis-s'l (4) 



where w{s) = dv/ds and where the integration is over a region for which A is zero at the 
boundaries. 

A widely used method involves replacing w with its asymptotic approximation This 
results in a convolution of the form, 

(,-t')4/3M^'), (5) 

or, 

_ 9/4 r ds' d\{s') 

"21737 (s_s')i/3 ds' ■ ^ ^ 

For this approach to be accurate, the characteristic length scale of the bunch must be 
much greater than This is satisfied for typical parameters of existing accelerators, 

but it is violated for highly compressed bunches, as well as for microbunched beams for 
proposed seeding schemes. An exact model for ID CSR is given in |5] and [6] that involves 
a more complicated wake function but that, in contrast with Eq. ([2]), does not assume small 
angles or small 7. In general, the computation of CSR convolutions is challenging for two 
reasons. First, the behavior of the wake function near s = becomes increasingly localized 
and large with increasing 7. Second, there can be a significant amount of numerical noise 
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FIG. 1. The single-particle wake, w{z), for a 100 MeV electron moving on a circle of radius 1 m, 
assuming large 7, small angles, and steady state. The full wake and its asymptotic approximation, 
w oc z~'^/^, are shown. The height of the full wake scales as 7^ and the width scales as 1/7^. 
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in the computed charge density, so methods that involve its derivative have to deal with 
this and usually involve careful smoothing. Below a new approach is presented that does 
not involve differentiating the charge density, and which furthermore includes short-range 
Wakefield effects. 

II. THE INTEGRATED GREEN FUNCTION APPROACH TO ID CSR 

A simple discretization of Eq. ([3| on a grid with cell size h leads to. 



where and Wk-k' denote the values of the line charge density and the wake function, 
respectively, defined on a grid of values Zk- This simple discretization has a potentially 
significant pitfall: It makes use of the wake function only at the grid points (even though it 
is known everywhere for certain problems where it can be determined analytically). This can 
be problematic when A and w have a disparate spatial variation, since, to achieve acceptable 
accuracy on a uniform grid, a sufficiently fine grid is needed to resolve the smallest features 
of both A and w. The use of a grid with variable grid spacing, or multiple grids of different 
resolution, can improve the efficiency of the algorithm, but with some additional algorithmic 
complexity. In the case of the CSR, the wake is extremely narrow and strongly peaked 
near z = 0. Depending on the bunch length and its internal structure (e.g., if it contains 
microbunches) , the asymptotic wake, the short-range wake, or both, might be important. 

There are numerous techniques for accurately and efficiently discretizing and evaluating 
Eq. ([3]). Qiang described a method for approximating convolutions to arbitrary accuracy 
using the Newton-Cotes formula [7j. For the CSR wake, Borland implemented a method 
of evaluating Eq. ^ that is used in the elegant code and other codes [8]. Carlsten im- 
plemented a method that includes some short-range radiation effects but requires solving 
a transcendental equation for the retarded times The approach below does not require 
solving a transcendental equation, and furthermore contains all short-range radiation effects 
contained in the single particle wake of Eq. ([T]). 

Integrated Green functions (IGF's) provide a means to accurately and efficiently compute 
convolutions when certain integrals involving the wake function can be obtained analytically 
or numerically [T014T2] . In the following we assume that the bunch length is so short com- 




(7) 
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pared to the bend radius that the longitudinal coordinate, 2;, can be used in place of the arc 
length, s. For our purposes we use linear basis functions to approximate the charge density 
within a cell of the computational domain. In that case Eq. ([s]) is approximated by. 



i,, Jo 



[{h - z')\k' + \k'+iz'] w{zk - Zk' - z'), 



(8) 



where the sum is over all values of k' in the discretization of A. Shifting the indices in the 
last sum above and collecting terms, the effective Green function is the coefficient of hX^' , 



k' 



(9) 



where 



w 



dz' [{h-z')w{C-z') + z'w{C + h-z')], (10) 

where C, = Zk — Zk' ■, and where we have used the notation w^-^-^ to denote w^^Jy. Rearranging 
terms we obtain. 



w 



eff 



/l2 



•C rC+h K K+h 

(h-C) I dz' wiz') + ih + C) / dz' wiz') + / dz' z'w{z') - / dz' z'w{z') 

C-h Jc JC-h Jc, 

(11) 



Integrating the last two terms by parts, and collecting terms, we obtain. 



w. 



eff 



where dv/dz = w. Finally, let 



dz' uiz) 



dz z/(z') 



C-h 



X{.z) 



dz'v{z'). 



It follows that 



w 



eff 
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(12) 

(13) 
(14) 



Note, however, that because integration by parts was used to arrive at this result, it will 
not necessarily hold for all values of (. In particular, for many problems w or dw/dz is 
discontinuous at z = 0. In such cases w^-^^ should be calculated separately. 

The integral in Eq. (13) can be done analytically for the ID GSR wake of Eqs. ([l])-([2]) 
(keeping in mind that the wake is zero for ^ < 0.). It follows that the effective Green 
function in Eq. ^ is given by. 
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eff 
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k,0 



{k > 0), 



(15) 



with w^J^ 



Xk 



for A; < 0. In the above, Xk is given for A; > by, Xk ■ 

[(1 - ^T)liA' 



§^Xk, where 
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+ log 



1 + + nth 



{k > 0), 



(16) 



where fifc = /x^ + ^y jj^l + 1 and /ifc = ^Zk- Also, Xfe = for < 0. Note that x is 
not singular, but, to evaluate it numerically for small argument one has the expansion, 
X=f;.{6-log(f) + ¥-l^---}for/i«l. 

The preceding specifies the effective Green function obtained from the full wake. If, 
instead of the full wake, its asymptotic expansion were used (as in Eq. (|5])), the effective 
Green function is still given by Eq. (15) but with x^ replaced by x^^™^' where 



" asymp 



27/8 

^V3 



f^T ik>0), 



(17) 



with Xk 



asymp 



OioT k < 0. 



III. EXAMPLES 



First consider a Gaussian line charge density with rms width a. The results that follow 
are proportional to the bunch charge, which we take to be 1 nC. The bend radius of curvature 
is p = 1 m. Figs. [2] and [3] show the calculated radiation electric field using the IGF approach 
obtained from the asymptotic wake and from the full wake. Fig. [2] has E = 100 MeV; the 
left, center, and right figures have a = 10, 1, and 0.1 micron, respectively. As is seen on 
the left, for the 10 micron case the asymptotic model gives results that are not significantly 
different from the exact model. The differences are larger for the 1 micron case in the center. 
In Fig. |2} right, the result based on the asymptotic wake has a large error. To understand 
this behavior, note that the characteristic length scale for this example is of order a, so the 
value of fi at the characteristic length scale is fichar ~ 37'^/(2p)a = 114., 11.4, and 1.14 when 
cr = 10, 1, and 0.1 micron, respectively. In the 0.1 micron case the condition fichar » 1 is 
badly violated, resulting in significant error. In contrast to Fig. [2] right. Fig. |3]also shows 
results with o" = 0.1 micron but at -E = 1 GeV. In this case, due to the higher energy 
fJ'char ~ 1100, so the asymptotic approximation is valid even for this short bunch length. 

The preceding demonstrated the ability of the IGF approach with the full wake function 
to capture the physics of short-range radiation effects. In addition, the IGF approach is more 
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FIG. 2. Radiation component of the azimuthal electric field based on effective Green functions 
derived from the asymptotic wake and from the full wake. The beam is a 1 nC, 100 MeV Gaussian 
bunch with rms width a. Left: cr = 10 micron. Center: a = 1 micron. Right: a = 0.1 micron. At 
a characteristic length scale of order a, fJ-char ~ 37^fT/(2p) = 114., 11.4, and 1.14 when a = 10, 1, 
and 0.1 micron, respectively. The condition fichar >> 1 is violated in the right-hand figure, leading 
to significant error in the GSR calculation based on the asymptotic model (red curve). 
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FIG. 3. Radiation component of the azimuthal electric field based on effective Green functions 
derived from the asymptotic wake and from the full wake. The beam is a 1 nG, 1 GeV Gaussian 
bunch with rms width a = 0.1 micron. 
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computationally efficient than the non-IGF approach when the density varies slowly com- 
pared with the wake function. Returning to the case with cr = 10 micron, Fig. |4] shows plots 
of the relative error in the radiation electric field, {E^^rad,num - E^^rad,exact) / E^^rad,exact,max, 
where Ez^rad,num IS the numerical result, Ez^rad,exact is the "exact" result based on the IGF 
approach with 500,000 grid points, and where the difference has been normalized by the 
maximum value of the "exact" field. Fig. |4| left shows the IGF, full wake results using 128, 
256, 512, 1024, 2048, and 4096 grid points on a domain covering ±10cr. The relative error 
is less than 0.1% with just 128 points. This good accuracy is explained by the fact that 
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128 grid points are sufficient to reasonably resolve the Gaussian charge density over this 
domain; it does not matter that 128 grid points are not sufficient to resolve the short-range 
structure of the wake. Fig. |4| right shows the non-IGF, full wake results using 1024, 2048, 
4096, 8192, and 16384 grid points. In this case 4096 grid points are required to achieve an 
accuracy of 0.1% or less. Taking a different viewpoint, at 1024 grid points the accuracy of 
the IGF approach is .0035% or less, while that of the non-IGF approach is 85% or less. 
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FIG. 4. Convergence test of the IGF method (left) and the non-IGF method (right) showing the 
relative error in the radiation component of the electric field for varying number of grid points. The 
beam is a 1 nC, 100 MeV Gaussian bunch with rms width cr = 10 micron. For these parameters 
the IGF method requires many fewer grid points (just 128) compared with the non-IGF method to 
achieve acceptable accuracy. Alternatively, for these parameters the IGF method is up to several 
orders of magnitude more accurate than the non-IGF method at a given number of grid points. 



Lastly, consider a Gaussian line charge density with sinusoidal microbunching, 

C 



X{z) 



-.-12. ^ ^ _ ^^^^ sin(fc^^))). 



(18) 



where C is a normalization constant, km = 27r/Am, Am is the microbunching wavelength, 
and where fmin controls the minimum depth of the modulation. Fig. |5] shows results with 
0" = 10 micron, A = 100 nm, and fmin = 0.1. The left figure is for an energy of 100 MeV, 
and the right figure is for 1 GeV. Both figures show IGF results with the full wake and with 
the asymptotic wake. Particularly noteworthy is the difference in microbunching-driven field 
enhancement between 100 MeV and 1 GeV: At 100 MeV the enhancement is approximately 
35%, at 1 GeV it is about 900%. Also, it is clear from the left plot that, for the 100 MeV 
parameters, the asymptotic IGF over-predicts the GSR field by a factor of 8, making the 
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results unusable; on the right at 1 GeV the asymptotic IGF gives results that are reasonably 
close to the full-wake results. This is consistent with the earlier observation, namely, that 
the asymptotic- wake IGF can be used for the 1 GeV parameters (since >> 1), but the 
full-wake IGF is required for the 100 MeV case. 
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FIG. 5. Radiation component of the azimuthal electric field in a Gaussian bunch with sinusoidal 
microbunching. The bunch rms width is cr = 10 micron, and the microbunching wavelength is 
Am = 100 nm. Left: 100 MeV. Right: 1 GeV. The red curve on the left is wrong because, for these 
parameters at 100 MeV, the asymptotic wake is not appropriate since fichar ~ 1- 

IV. DISCUSSION 



A new convolution-based method has been presented for computing the ID GSR radia- 
tion electric field. Using an integrated Green function (IGF) approach, an effective Green 
function was derived based on the ID GSR wake. If the full single-particle wake of Eq. ([T]) is 



used, the resulting effective Green function is given by Eqs. (15) and (16). If the asymptotic 
wake is used, the resulting effective Green function is given by Eqs. ( [I5| and (17). In con- 
trast with methods that use the derivative of the line charge density dX/dz, in this approach 
the discrete convolution is performed using A. If the analysis starts from the full GSR wake, 
the resulting convolution is able to capture both short- and long-range effects. A numeri- 
cal example demonstrated the ability of the full-wake IGF approach to model short-range 
effects when the asymptotic model is not appropriate. The example also showed that the 
IGF approach has the potential to achieve accurate results, but with far fewer mesh points, 
N, than the non-IGF approach, when A varies on a scale that is long compared with that of 
the single-particle wake. The ability to use fewer mesh points has two consequences: first, 



an IGF code will run faster since the execution time scales as 2A^log(2A^) for an FFT-based 
method, and second, there will be less numerical noise since there will be more particles per 
cell. (Note, however, that if the scale of the bunch microstructure less than or comparable to 
that of the short-range wake, then the mesh points needed to resolve A will also resolve the 
short-range wake, and there might be no advantage to using an IGF.) Lastly, an example 
with microbunching showed significant enhancement in a IGeV bunch with microbunching 
at a wavelength of 100 nm. This enhancement might (as in Fig. [sj right) or might not (as 
in Fig. [sj left) be modeled correctly in a code based on the asymptotic convolution, Eq. ^ 
or Eq. ([6]), depending on the problem parameters. The validity of the asymptotic model 
is most likely to be of concern at low energy, and in short bunches or in bunches with a 
microstructure. 

The accuracy of the IGF approach depends on choosing a grid that is sufficient to resolve 
the charge density, not the single particle wake function. Assuming that is the case, an IGF- 
based particle-in-cell code will usually produce acceptable results. In contrast, a non-IGF 
code can produce unusable results, even when the grid is fine enough to resolve the charge 
density, if the Green function is not well-resolved by the grid. In general IGF codes tend to 
be more robust than their non-IGF counterparts since they are less sensitive to resolution 
issues associated with the grid spacing. 

The approach described here does not require computing the derivative of the line charge 
density, dX/dz. This is desirable given the loss of accuracy in computing the numerical 
derivative of a possibly noisy, numerically derived function. In principle one could make use 
of d\/dz in an IGF code by applying the IGF technique to Eq. Q. But there is no obvious 
advantage to doing this since, whether one uses Eq. ([3]) or (|4]), the singularity is handled 
analytically. From a programming standpoint there is never any need to compute dXjdz. 
Instead one simply deposits charge on a ID grid and performs the discrete convolution given 



by Eqs. (15) and (16). In regard to computational effort, numerical convolutions can be 
computed rapidly by turning the discrete convolution into a cyclic discrete convolution that 
can be handled using Fast Fourier Transforms methods (see the Appendix of jl3j). 

The examples presented illustrate the dependence of the radiation field on both energy 
and bunch length scale: On one hand, the short-range wake grows as 7^, so there is potential 
for huge field enhancement at high energy, with important implications for microbunched 
beams. On the other hand, if ~ T'sj p » 1 then the short-range portion of the wake is 
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irrelevant, the asymptotic portion dominates, and the radiation field is independent of 7. 

We have presented an IGF-based method to calculate ID GSR. In light of that it is 
reasonable to ask, under what conditions is the ID approximation likely to be valid? An 
approximate criterion can be obtained as follows: The condition that the radiation cones of 
two electrons separated by a transverse distance x overlap is, for small angles, x ~ pOipret 
where p is the bend radius, 6 is the opening angle, and where ipret is the retarded angular 
position. For frequencies smaller than the critical frequency Uc = (3/2)7^c//9, the coherent 
radiation is primarily within the angle 9 = {1/'j)(uJc/ujY^^- The dominant radiation is at 
frequencies of order of the inverse bunch length, u ~ 27rc/cr. The range of retarded angles is 
related to the bunch length by ijjret ~ {2Aa/ pY^^. It follows that the overlap condition for a 
bunch of rms transverse size x and rms length o" is x ~ p{a/pY^^, and the ID approximation 
is likely to be valid as long as x is well below this value. We have found this criterion to be 
consistent with results from a particle-based, massively parallel Lienard-Wiechert code |T3] . 

Lastly, two phenomena should be mentioned that are not included in this analysis. First, 
we have discussed just the radiation component of the field, not the Coulomb component. 
Since the Coulomb field varies as l/r^, it can't be handled in the one-dimensional approxi- 
mation since the integral in Eq. (|3| would be divergent. It can, however, be handled using a 
two-dimensional model. Second, we have not discussed shot noise effects. Simulations using 
a massively parallel Lienard-Wiechert code have shown that shot noise effects can be a major 
consideration near 1 GeV and above, even in bunches with several billion electrons [H]. In 
general, classical CSR effects in a bunch containing a finite number of particles involve a 
smooth component and a stochastic component. The analysis presented here describes an 
efficient way to model the smooth component of the radiation field in ID, and to include 
both short- and long-range phenomena if needed. 



In conclusion, the method described by Eqs. (15) and (16) provides a robust method 
for computing ID CSR effects. Since it includes both short- and long-range effects, it will 
produce more reliable results than the asymptotic model after the microbunching instability 
sets in, {i.e., after small-scale structures develop), resulting in simulations that provide 
more reliable estimates of the microbunching growth rate. Given the simplicity of the 
algorithm, it's numerical efficiency, and the added benefit that it does not require numerically 
differentiating the line charge density, this provides an attractive new method for modeling 
ID steady-state CSR effects. 
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